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LONG-TERM  GOALS 

Our  proposed  research  will  focus  on  the  horizontal  underwater  acoustic  communications  (UAC).  One 
important  metric  for  evaluating  UAC  systems  is  the  maximum  achievable  range-rate  product.  Building 
on  co-located  and/or  distributed  multiple  transducers  and  hydrophones,  we  plan  to  develop  rate- 
oriented  and  range-oriented  multiple-input  multiple-output  (MIMO-)UAC  systems,  as  well  as 
additional  range/rate  enhancement  via  the  deployment  of  distributed  transceiver  units  which  operate  in 
a  cooperative  manner. 

OBJECTIVES 

Coherent  underwater  acoustic  communication  is  of  great  importance  for  high  data  rate  applications, 
such  as,  image  transmission  or  multiuser  networks  at  long  ranges.  We  consider  MIMO  transmit 
beamforming  under  the  uniform  elemental  power  constraint.  This  is  a  non-convex  optimization 
problem,  and  it  is  usually  difficult  to  find  the  optimal  transmit  beamformer.  We  will  first  find  a 
beamforming  solution  in  an  ideal  setup,  and  then  consider  practical  finite-rate  feedback  methods 
required  to  implement  the  transmit  beamforming. 

The  time-varying  and  frequency-selective  underwater  acoustic  channel  consists  of  a  great  challenge  for 
channel  estimators  necessary  for  coherent  communications.  Such  difficulty  is  aggravated  when 
multiple  transducers  and/or  hydrophones  are  deployed.  Therefore,  we  also  investigate  differential 
MIMO  transmission  and  reception  approaches  tailored  for  doubly-selective  MIMO  channels.  To 
enhance  the  system  error  performance,  we  notice  that  the  three-dimensional  space-time-frequency 
variation  of  the  underwater  channel  can  be  exploited  to  provide  three-dimensional  diversity  gain. 

Based  on  all  these,  our  second  objective  is  to  develop  a  differential  MIMO  transceiver  tailored  for 
doubly-selective  underwater  channels. 

APPROACH 

1.  MIMO  Transmit  Beamforming  under  Uniform  Elemental  Power  Constraint 

The  Cyclic  Algorithm  (or  MIMO  Transmij  Beamformer  Design 

Consider  an  (NnNr)  MIMO  communication  system  in  a  quasi-static  frequency-flat  fading  channel.  At 

r  ir 

the  transmitter,  the  symbol  5  e  C  is  modulated  by  the  beamformer  w,  =  wu  w, ,  ...  w;  v  !  .  At  the 
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receiver,  after  processing  with  the  combining  vector  wr  =  j  wr ,  wr  2  . . .  wr  N  ,  the  sampled 

baseband  signal  is  given  by  y  =  w'  [Hw,s  +  n],  where  HeCw  is  the  channel  matrix.  Our  problem 
at  hand  is 

max  |w;Hw,  |2,  subject  to  |  w,(.  |2=-!-,  i  =  \,2,..^Nt,  ||wr||  =  l.  (1) 

•  yv( 

This  is  a  non-convex  optimization  problem,  which  is  usually  difficult  to  solve,  and  no  globally  optimal 
solution  is  guaranteed.  We  solve  it  in  the  following  way: 

First,  we  consider  the  multi-input  single  output  (MISO)  case  and  propose  a  cyclic  algorithm: 

•  Step  0:  Set  wr  to  an  initial  value  (e.g.,  the  left  singular  vector  of  H  for  its  largest  singular  value). 


•  Step  1:  Obtain  the  w(  that  maximizes  (1)  for  wr  fixed  at  its  most  recent  value.  By  taking  w’H  as 


the  “effective  MISO  channel,”  the  optimal  solution  is:  w,  =  . —  e 

yjN, 


jz  hV 


•  Step  2:  Determine  the  wr  that  maximizes  (1)  for  w,  fixed  at  its  most  recent  value.  The  optimal  wr 


is  the  MRC:  wr 


Hw, 

|Hw, 


.  Iterate  Steps  1  and  2  until  a  given  stop  criterion  is  satisfied. 


Finite-Rate  Feedback  (or  TransmU  Beamformins  Designs 

We  consider  herein  vector  quantization  under  uniform  elemental  power  constraint  (VQ-UEP)  for 
finite-rate  feedback.  In  this  case,  both  the  transmitter  and  the  receiver  maintain  a  codebook  which  is 
used  to  reconstruct  the  transmit  beamformer.  The  generalized  Lloyd  algorithm  is  usually  implemented 
for  codebook  construction.  The  design  criterion  is  adopted  to  satisfy  two  conditions:  the  nearest 
neighborhood  condition  (NNC)  and  the  centroid  condition  (CC).  NNC  is  to  find  the  optimal  partition 
region  for  a  fixed  codeword,  while  CC  updates  the  optimal  codeword  for  a  fixed  partition  region.  The 
monotonically  convergent  property  is  guaranteed  due  to  obtaining  an  optimal  solution  for  each 
condition.  We  construct  the  codebook  offline  by  maximizing  the  average  receive  SNR: 

First,  we  generate  a  training  set  from  a  sufficiently  large  number  of  channel 


realizations.  Next,  starting  from  an  initial  codebook,  we  iteratively  update  the  codebook  according  to 
the  following  two  criteria  until  no  further  improvement  is  observed. 

•  NNC:  for  given  {w,}j!v, ,  assign  H,;  to  the  i  th  region  S.  =  {H„ :  ||H„wJ2  >  ||H„wyf  ,Y/  *  /},  where 
S,,i  =  1,2,..,  Nv,  is  the  partition  set  for  the  i  th  codeword  w,  • 


CC:  for  given  partition  St ,  the  updated  {w,},^  satisfy  w^argmax/T,  ||Hn 

to  I H'l.m |2  =  ~r »  /M  =  l,..,Af„for  /  =  1,2, Let  R,- =£<I[H’H  |H  e 

N, 


w|2 1  H„  eSj, subject 
S',]  and  R,1 2  be 


Hermitian  square  root  of  R;. 

For  a  given  codebook  W  :=  {wdWjv »Watv}>  the  receiver  first  chooses  the  optimal  beamformer  as: 
Wopt  =  argmax  II Hvvll  .Then  we  feedback  the  index  of  Vvopt  to  the  transmitter  with  log,  N  -  B  bits. 
Average  Degradation  of  the  Receive  SNR 

We  quantify  below  the  effect  of  VQ-UEP  with  finite-bit  feedback: 

Proposition  1  For  i.i.d.  MISO  Raleigh  fading  channels ,  the  average  degradation  of  the  receive  SNR, 
for  an  Nt  -antenna  transmit  beamforming  system  with  an  Nv  =  2s  -size  VQ-UEP  codebook,  can  be 
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approximated  as: DV(B)  ~  (Nt  —  1)  •  2l 


(2~b  +  (1  -  a)"'-1)"'"1  -  (1  -  a)Nt 


al  -Nt(  1  -  a)o£. 


and  ofi  is  the  covariance  of  each  scalar  channel. 


,  1  ,  7T  Nt 

where  a  =  — - 

Nt  4(Wt-l) 

The  average  degradation  of  the  receive  SNR  can  be  proven  to  be  monotonically  decreasing  with 
respect  to  non-negative  real  number  B .  Given  a  d  egradation  amount  D0 ,  this  proposition  provides  a 
guideline  to  determine  the  necessary  number  of  feedback  bits. 


2.  Orthogonal  Space-Time  Block-Differential  Modulation  over  Underwater  Acoustic  Channels 

Our  approach  hinges  upon  the  delay-Doppler  spread  channel  characterization  and  a  block-by-block 
transmission  and  receiver  processing.  Without  loss  of  generality,  let  us  consider  a  SISO  doubly- 
selective  UAC  channel  with  impulse  response  h(t;  r)  and  frequency  response  H(f;  t)  which  is  also 
known  as  the  delay-Doppler  spread  function  [1],  The  I/O  relationship  of  such  a  channel  can  then  be 

obtained  as  y(t)  =  flmax  x(t  —  r )H(f;  r)e,27r^t  drdf,  where  Tmax  is  the  maximum  delay 

spread  and  fmax  is  the  maximum  Doppler  spread.  Consider  a  block-by-block  transmission  with  N 
symbols  per  block  and  symbol  duration  Ts.  It  follows  that  the  time-  and  frequency-domain  resolutions 
of  the  corresponding  discrete-time  equivalent  system  are  Ts  and  (NTs)-1,  respectively.  Hence,  the 

.  2  nk  . 

discrete-time  equivalence  of  the  above  I/O  relationship  is:  y(n)  =  Y.L=oY,k=-K  x(n  ~  0  H(k;  l)eJ  »  , 

L  =  [“"p"] > K  =  \fmax  NTS 1-  Denoting  Q  =  2 K,  to  ■=  2n/N  and  hq(l)  =  H(q  -  K;  l),  we  have  [2]: 

y(n)  =  Y.i=ox(n  ~  0  In  short,  we  obtain  a  basis-expansion  model  (BEM)  as: 

Q 

h(n;  0=2>,  (l)e>  ^-K)a)n  ,l€[0,L].  (3) 

<j=0 

Then  the  system  I/O  relationship  becomes  y(n)  =  £f=o  h(n;  l)x(n  —  l)  +  z(n),  or  y  =  Hx  +  z. 

For  a  frequency-selective  but  time-invariant  channel,  H  is  Toeplitz,  which  can  be  diagonalized  by  FFT 
and  IFFT  operations  together  with  cyclic  prefix  (CP)  to  enable  differential  encoding.  In  the  time- 
varying  case,  however,  the  elements  on  (sub)diagonals  are  different.  Fortunately,  based  on  the  BEM 
model,  one  could  diagonalize  the  channel  matrix  via  two  FFT-like  operations,  one  deals  with  the 
multipath-frequency  conversion  as  in  conventional  OFDM,  the  other  the  time-Doppler  conversion 
(n  —  q  in  (3)).  This  diagonalization  will  then  enable  differential  coding.  Special  considerations  will  be 
added  to  enable  Doppler-multipath  diversity.  On  top  of  these,  we  also  employ  the  differential 
orthogonal  space-time  block  code  (OSTBC)  [3],  [4]  to  harvest  the  space  diversity. 

WORK  COMPLETED 

1.  MIMO  Transmit  Beamforming  under  Uniform  Elemental  Power  Constraint 

First,  we  consider  the  bit-error-rate  (BER)  performance  of  our  beamformer  with  perfect  CSI  at  the 
transmitter.  For  comparison  purposes,  we  also  implement  the  conventional  transmit  beamforming 
design  without  the  uniform  elemental  power  constraint  (Con  TxBm)  and  the  conventional  design  with 
peak  power  clipping  (TxBm  with  Clipping),  in  addition  to  our  own  (UEP  TxBm).  Fig.  1  shows  the 
BER  comparison  for  (4,2)  and  (8,8)  MIMO  systems.  Con  TxBm  achieves  the  best  performance  since 
it  is  not  under  the  uniform  elemental  power  constraint.  Under  the  uniform  elemental  power  constraint, 
UEP  TxBm  schemes  outperform  TxBm  with  Clipping.  We  consider  next  the  effect  of  vector 
quantization.  We  use  m-bit  VQ-UEP  with  Nv  =  2m  ,  m=2,4,6,8  to  quantize  the  transmit  beamformer 
for  a  (4,  2)  MIMO  system.  Fig.  2  shows  that  VQ-UEP  outperforms  TxBm  with  Clipping  when  B  =  4,  6, 
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8,  although  the  latter  is  with  the  perfect  feedback.  By  using  the  6-bit  VQ-UEP,  the  (4,2)  MIMO  system 
only  suffers  from  slight  performance  loss  compared  to  UEP  TxBm  with  perfect  feedback.  When  more 
bits  are  used,  we  can  further  close  the  small  gap.  In  Fig.  3,  we  carry  out  Monte-Carlo  simulations  for  a 
(4,1)  system  and  plot  the  numerically  simulated  degradation  result.  We  can  see  that  the  closed-form 
expression  in  Proposition  1  gives  a  very  accurate  approximation. 

2.  Orthogonal  Space-Time  Block-Differential  Modulation  over  Underwater  Acoustic  Channels 

Eq.  (3)  implies  that  the  N(L  +  1)  time-varying  channel  parameters  h(n>  l)  can  be  captured  by  the 
(Q  +  1  )(L  +  1)  time-invariant  (within  one  block)  channel  parameters  hq  (/).  In  practice,  however,  we 
could  only  observe  h(n\  /).  The  first  problem  at  hand  is  how  to  form  hq(l)  from  the  observed  h(n;  l) s. 

Lemma  1  With  Wb  =  [a)b, a)b+B,  ...,a)b+QB]  and  N  =  B  x  (Q  +  1),  we  have:  WblVb  =  (Q  +  1) 

/<?+ 1;  i.  9(0  =  ^~1  <  [h(b;  o,  Kb  +  B;  I) . h(b  +  QB ;  l)]T  where  g(l)  =  [hod). ...  hQ(Q)]T. 

Notice  that  Wb  is  reminiscent  of  the  FFT  matrix.  This  lays  the  foundation  for  our  2-dimensional  FFT- 
based  differential  approach,  which  is  detailed  in  the  following  steps. 

Step  7;  CP  Insertion  and_  Removal 

At  the  transmitter,  every  frame  is  partitioned  into  P  subblocks  each  containing  ( Q  +  1  )M  symbols. 
Every  M  symbol  segment  is  augmented  by  a  cyclic  prefix  (CP)  of  length  L,  which  is  removed  at  the 
receiver.  The  transmitted  symbol  frame  structure  is  shown  in  Fig.  4.  As  a  result  of  this  step,  the 
channel  matrix  H  becomes  block-diagonal  (see  Fig.  5). 

Step  2:  Interleaving.  and_  De-Interleavins 

Lemma  1  implies  that  recovery  of  hq(l)  using  Wb  requires  equally-spaced  h(n,  l) s.  Hence,  we 
perform  interleaving  and  de-interleaving  at  the  transmitter  and  receiver,  respectively:  x  :=  07  x  and 
y  :=  0y ,  where  Q/Q'  takes  every  PM  columns/rows.  As  a  result,  the  I/O  relationship  becomes 
y  =  Hx  +  z,  where  H  :=  QHQT  consists  of  block-quasi-circulant  matrices  on  its  diagonal  as  depicted 
in  Fig.  6.  Because  of  the  diagonal  structure  of  H ,  the  I/O  relationship  can  be  represented  by  shorter 
streams:  y(p)  =  W(p)X(p)  +  z(p),  p  E  [0,  P  -  1], 

Step  3:  Time-Doppler  FFT 

Now  we  have  an  equivalent  channel  matrix  H(p )  which  consists  of  diag {h(n>  0}-  Lemma  1  indicates 
that,  in  order  to  obtain  g(l ),  we  need  to  multiply  Wb  with  h(n ;  /);  that  is,  to  multiply 

Wp(M+L)+L  °0?+i)x(<?+i)  -  0((?+1)x(<?+1) 

Sl(p)  :=  °«?+l)x(<?+l)  ^P(M+L)+Z,+1 

.0«?+l)x(<?+l)  0«?+l)x«?+l)  ^p(M+i)+t+M- 1. 

with  y(p)  at  the  receiver.  In  the  meantime,  x(p)  =  U(p)  <8>  l<?+i  is  also  needed  at  the  transmitter, 
where  t/(p)  is  the  M  x  1  data  vector.  Then,  the  I/O  relationship  becomes  y(p)  =  (Q  +  1  )GU(p)  + 
z(p),  where  G  is  the  M(Q  +  1)  x  M  equivalent  channel  matrix  (see  Fig.  7),  which  is  independent  of  p. 

Step  4:  Multipath-Frequency  FFT 

It  is  well  known  that  IFFT  at  the  transmitter  and  FFT  at  the  receiver  can  diagonalize  a  circulant  matrix. 
The  transmitted  vector  can  be  formed  as  l/(p)  =  Fj^u(p),  where  FM  is  the  M-point  FFT  matrix.  And 
the  received  vector  v(p)  :=  FMy(p),  with  FM  ~  [lQ+1  ®  fx,  ...,IQ+1  ®  fM\.  Therefore  v(p)  = 

Du(p )  +  t]{p)  is  now  the  I/O  relationship  with  D  being  a  block-diagonal  matrix  containing  the 
frequency  responses  of  channel  vectors  g^  ■=  [h9  (0), ... ,  hq  (L)]  ,  as  depicted  in  Fig.  7. 
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Step  5j  MIMO  Differential  Encoding 

It  is  now  straight  forward  to  do  differential  encoding  if  the  system  is  single-input  single-output  (SISO). 
When  multiple  transducers  are  deployed,  we  extended  the  differential  OSTBC  [3],  [4]  for  time- 
invariant  flat-fading  channels  to  doubly-selective  channels.  Without  loss  of  generality,  we  considered 
the  case  consisting  of  Nt  =  2  transducers  and  Nr  =  1  hydrophones.  We  group  2  consecutive 
subblocks  into  one  and  index  it  with  pb.  The  differential  encoding  goes  as:  [ui(Pfc),“2(Pfc)]  = 

rui(pi-l)  -u2(P/,-i)* 


[Ddl  (Pb).  f>d2(Pi)]  [ 


«2(Pft-l)  «l(Pb-l)’ 


,  where  Ddl(pb )  and  Dd2(pb )  are  diagonal  data  matrices, 

Ui  ( pb )  and  u2(pb)  are  from  the  p^-th  block.  These  are  then  mapped  to  different  transducers 
(subscript)  as:  =  Ui(Pi>).Ui(2pb  +  1)  =  ~n2(_pby and  u2(2pb)  =  u2(pi),u1(2pi  +  1)  = 

Ui(pt)‘.  With  v(2pb)s  denoting  the  received  blocks,  the  ML  decoder  is 

v(2pbM2pb  -  iy  +  v(2pb  +  l)*v(2(p„  -  1)  +  1) 

--Om(pft) 


0ci(pi,)  =  argmin 
Do  l 


K2(pfc-l))||  +  \\v(2(pb  -  1)  +  1)||2 


Simulations  and_  comparisons: 

Putting  everything  together,  the  transceiver  structure  is  shown  in  Fig.  8.  Using  this  structure,  we 
performed  simulations  with  three  types  of  channels:  time-invariant  (TI),  time-varying  (TV)  and 
experiment-extracted  (EXP)  (see  Fig.  10).  Fig.  9  depicts  the  diversity  differences  among  the  BER 
curves.  Comparing  the  two  curves  using  QPSK  and  k  =  3  coding  over  TI  channels,  we  see  that  the 
coding  in  [5]  provides  diversity  advantage.  Doppler  diversity  can  be  observed  by  comparing  the  TI  and 
TV  cases.  From  the  curve  generated  using  the  SISO  scheme  in  [6]  and  the  one  using  our  scheme  over 
TV  channels,  the  spatial  diversity  advantage  is  evident.  Fig.  9  also  shows  the  performance  degradation 
encountered  by  plain  OFDM.  With  experiment-extracted  channels,  the  BER  curve  exhibits  an  error 
floor  due  to  the  noise  in  the  estimated  channels.  Recall  that  the  BEM  has  a  requirement  on  fmax  and 
the  noise  spikes  introduce  high  frequency  components.  Therefore,  this  heavy  tail  is  expected  to 
disappear  in  real  experiments. 


RESULTS 


Our  MIMO  transmit  beamformer  accounts  for  the  uniform  elemental  power  constraint,  which  is  critical 
for  the  power  efficiency  and  lifetime  of  high-rate  UAC  transmitters.  While  this  problem  is  a  difficult 
non-convex  optimization  problem,  we  managed  to  develop  a  high-performance  low-complexity  cyclic 
algorithm.  We  also  considered  practical  finite-rate  feedback  and  studied  the  average  SNR  degradation. 
Although  the  channel  assumptions  we  employed  still  differ  from  realistic  underwater  acoustic  channels, 
our  methods  can  be  readily  extended  to  the  underwater  acoustic  applications. 

Doubly-selective  MIMO  UAC  channels  are  ample  in  diversity  but  are  difficult  to  estimate.  We  are 
among  the  first  to  explore  full-diversity  differential  approaches  in  such  channels.  Our  simulation 
results  using  statistically  modeled  and  experiment-extracted  channels  are  promising.  In  addition  to 
field  tests,  we  will  also  continue  working  on  improve  the  data  rate  of  our  diversity  differential 
transceivers. 


IMPACT/APPLICATIONS 

The  natural  bandwidth  limitations  of  coherent  underwater  acoustic  channel  suggest  a  technical 
breakthrough.  MIMO  signal  processing  is  a  promising  bandwidth  efficient  method  to  high  data  rate 
and  high  quality  services.  Taking  both  coherent  and  differential  approaches,  our  promising  results  are 
expected  to  favorably  impact  high-rate  long-range  MIMO-UAC  designs. 
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Figure  1:  Performance  comparison  of  various  transmit  beamformer  designs  with  perfect  CSI  at  the 

transmitter:  (a)  the  (4,2)  case,  and  (b)  the  (8,8)  case. 

[graph:  Our  design  (red)  is  much  better  than  existing  approach  (blue),  and  is  dose  to  the  design 

without  uniform  power  constraint  (black).] 
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Figure  2:  Performance  of  vector  quantization  under  uniform  elemental  power  constraint  with 
various  bits  feedback,  for  a  (4,2)  system.  Note  that  8-bit  VQ-UEP  and  UEP  TxBm  with  perfect 
feedback  curves  almost  coincide  with  each  other. 

I graph :  With  8-bit  for  channel  feedback,  little  performance  degradation  observed.! 


Figure  3:  Average  degradation  of  the  receive  SNR  for  a  (4,1)  M1SO  system. 

I graph :  Our  theoretical  prediction  of  feedback  SNR  loss  agrees  closely  with  simulations.! 
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Figure  4.  Transmitted  symbol  frame  structure 
[graph;  A  transmitted  frame  with  P  =  2,  Q  =  2] 


Figure  5.  Sub-transmitting  matrix  structure 
[graph:  quasi-circulant  structure  of  each  sub-channel  matrix,  M  x  M/ 


Figure  6.  I/O  after  (de-)interleaving 

I graph :  quasi-circulant  structure  of  each  sub-channel  matrix,  M(Q  +  1)  x  M(Q  +  1)/ 
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Figure  7.  Multipath-Frequency  FFT 
[graph:  equivalent  I/O ,  block-circular ,  circular ,  block-diagonal] 
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Figure  8.  Overall  transceiver  structure 

[graph:  transmitter:  differential  encoding,  serial  to  parallel,  LCF  encoding,  I  FFT,  interleaving, 

insert  CP; 

receiver:  remove  CP,  de-interleaving,  FFT,  differential  decoding  / 
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Figure  9.  Bit  error  ratio  vs.  signal  to  noise  ratio 
[graph:  our  scheme  has  the  best  performance,  achieving  3-dimensional  diversity,  traditional 
OFDM  degrades  over  time-varying  channels,  our  scheme  also  perform  well  over  real  channels] 
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Figure  10.  One  snapshot  of  the  underwater  acoustic  channel  extracted  from  the  experiment 
I  graph:  Amplitude  vs.  delay  index,  the  channel  delay  may  as  long  as  40  taps] 

(Courtesy  Dr.  Daniel  Kilfoyle) 
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